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Transition voltage spectroscopy (TVS) has recently been introduced as a spectroscopic tool for 
molecular junctions where it offers the possibihty to probe molecular level energies at relatively low 
bias voltages. In this work we perform extensive ab-initio calculations of the non-linear current 
voltage relations for a broad class of single-molecule transport junctions in order to assess the 
applicability and limitations of TVS. We find, that in order to fully utilize TVS as a quantitative 
spectroscopic tool, it is important to consider asymmetries in the coupling of the molecule to the 
two electrodes. When this is taken properly into account, the relation between the transition voltage 
and the energy of the molecular orbital closest to the Fermi level closely follows the trend expected 
from a simple, analytical model. 
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Molecular electronics holds the promise of continu- 
ing the miniaturization of electronic devices beyond the 
limits of standard silicon technologies by using single 
molecules as the active elements [l|, l2[. In the design 
and characterization of molecular devices, the electronic 
structure of the molecule is naturally of vital impor- 
tance. In fact, when quantum inter fererence effects are 
disregarded [3j, there is a direct relation between a junc- 
tion's transport properties and the distance of the molec- 
ular energy levels to the electrode Fermi level. The latter 
may in principle be determined from peaks in the dJ/dF 
curve, where / is the current and V is the bias voltage. 
Assuming that the molecular level closest to the Fermi 
level is the highest occupied molecular orbital (HOMO) it 
would require a voltage oiV ^ 2|£^i? — Sh\ to probe the 
HOMO position in a symmetric junction. However, in 
practice the molecular junction often becomes unstable 
and breaks down due to the large current density before 
the peak in the d//dV is reached. 

Transition voltage spectroscopy (TVS) was introduced 
as a spectroscopic tool in molecular electronics by Beebe 
et al. y . Beebe found that the Fowler-Nordheim graph 
of a molecular junction, i.e. a plot of ln(//F^) against 
1/y, showed a characteristic minimum at a bias volt- 
age l^in which scaled linearly with the HOMO energy 
obtained from ultraviolet photo spectroscopy (UPS). Im- 
portantly, the TVS minimum is obtained at relatively 
low bias voltage before electrical break-down. TVS is 
now becoming an increasingly popular tool in molecular 
electronics 1:2, ignSj. 

The original interpretation of TVS introduced by 
Beebe [^ and applied in most later experimental work, 
is based on a Simmons tunnel barrier model [10| . Within 
this interpretation, the TVS minimum is obtained, when 
the tunnel barrier, due to the applied bias potential, 
changes from being trapezoidal to triangular. The tran- 



sition voltage equals the barrier height, which is inter- 
preted as the distance from the Fermi level to the closest 
molecular level. However, it has recently been pointed 
out by Huisman et al. that the barrier model is incon- 
sistent with experimental data jn| which on the other 
hand is more appropriately described by transport via 
a single electronic level. In the single-level model the 
transmission function was assumed to have a Lorentzian 
shaped, and all non-linear effects due to the finite bias, 
were neglected. While these assumptions may be reason- 
able, they are not obviously fulfilled in a realistic molec- 
ular junction at high electric fields. It is therefore of 
interest to compare the simple Lorentzian transmission 
model with more realistic calculations. 

In this Letter, we present a quantitative analysis of 
TVS based on extensive ab-initio calculations of the non- 
linear current- voltage relations for a broad class of molec- 
ular junctions. The ratio of the TVS minimum to the 
HOMO level position is found to vary between 0.8 and 
2.0 depending on the junction asymmetry, i.e. quite dif- 
ferent from the one-to-one relation assumed so far. The 
large variation is due to the difference in the non-linear 
response of the molecular level to the bias voltage. The 
importance of asymmetry effects is further signified by 
the fact that many of the experiments using TVS were 
performed on asymmetric molecules in an asymmetric 
conductive AFM measurement setup [j-lsl, lS|- Indeed, 
as we show below a larger degree of consistency between 
TVS measurements and UPS data is obtained when the 
asymmetry is taken into account. 

We begin our analysis by considering a simple model 
for transport via a single electronic level equivalent to the 
work in Ref. [ll|. The transmission function is assumed 



to be a Lorentzian 
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Figure 1: Ratio between molecular level energy and transition 
voltage, 1^0 — EF\/Vmin, vs. molecular level energy, so and 
broadening, F, for a symmetric junction, ?7 = 0, (left), and a 
completely asymmetric junction, 77 = 1/2, (right). Note the 
different scales of \£o — EF\/Vmin for the symmetric and asym- 
metric junctions. The lower part illustrates how the molecu- 
lar level moves under a finite bias voltage. At zero bias the 
level is located at £0, but when a bias voltage is applied, the 
level follows to some degree the chemical potential of the left 
electrode, due to a stronger coupling to this electrode. 



of the current depends quadratically on the bias voltage, 
/ (X V^, and d{\n{I/V'^))/dV = 0. An example of a cal- 
culated Fowler- Nor dheim plot can be seen in the inset of 
Fig. m In Fig. [T]we plot the ratio l^o — ^F|/^min as a 
function of broadening, F and molecular level energy, Sq, 
for a symmetric junction (77 = 0) and a completely asym- 
metric junction {r] = 1/2). Fig. [1] illustrates two main 
points: First, the ratio l^o — ^^Fl/'Knin is nearly constant 
over a large range of F and sq values. For a given de- 
gree of asymmetry (value of 77), the TVS minimum can 
therefore be used as a direct measure of the molecular 
level position, independently of F. Note, however, that 
for F/I^o — Ef\ > 1, there is no minimum in the Fowler- 
Nordheim plot, since the molecular level is too close to 
the Fermi level, and the current increases slower than 
oc V'^ at all bias values. The second conclusion from Fig. 
[T]is that the ratio \eo — Ef\ /Vmin ranges from 0.86 to 2.0 
depending on the asymmetry of the molecular junction. 
In order to use TVS as a quantitative tool, knowledge of 
the asymmetry factor is therefore needed. On the other 
hand, if the molecular levels can be determined by other 
means, the TVS can by used to measure the asymmetry 
factor. 



where eo and F are the molecular level energy and broad- 
ening, respectively. / is a constant factor to account for 
multiple molecules in the junction, asymmetric coupling 
to the electrodes, etc. Our analysis will not be depen- 
dent on the actual value of /. We shall assume that the 
molecular level is the HOMO level, i.e. Eq < Ep. This 
is the typical case for thiol bonded molecules (as consid- 
ered in this work) , where electron transfer from the metal 
to the sulphur end group shifts the molecular levels up- 
ward in energy [l2|. At finite bias voltage, the molecular 
level may be shifted relative to the zero-bias position, as 
shown schematically in Fig. [1] (lower part). Such non- 
linear effects express themselves differently in symmetric 
and asymmetric junctions [l3|, [l^. For example, if the 
molecule is only strongly coupled to the left electrode, 
the molecular levels will follow the chemical potential of 
the left contact. On the other hand, for a symmetric 
junction the molecular level will remain at the zero-bias 
position. We describe the degree of asymmetry by the 
parameter r] G [—1/2; 1/2] such that the current is given 

by 
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T{E;eo^vV,r)[h{V)-fR{V)]dE, (2) 



where fL/niV) = l/[exp {Ef ± eV/2)/kBT + 1] are the 
Fermi-Dirac distributions for the left and right contact, 
respectively. The dependence of the level position on the 
bias voltage is the main difference between our model and 
the one considered in Ref. lul. 

Using Eq. (|2]) we calculate ln(//'\/^) and find the min- 
imum at the bias voltage Fmin- At this voltage, the slope 
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Figure 2: Ratio between the HOMO energy (at zero bias) 
and the transition voltage, Vmin, vs. asymmetry parameter, r\. 
Solid line is obtained from a Lorentzian transmission function 
using Eqs. ([T]) and ([2|) and symbols are results of ab-initio 
finite bias calculations. 

We now turn to our ab-initio finite bias calculations 
and Fig. [21 which is our main result. It shows the ratio 
of |6:h — EF\/Vmin VS. T] for the 17 molecular junctions 
listed in Fig. [3l The solid line is the result obtained 
from the one- level model using Eqs. ([T]) and (J2j). Al- 
though there are deviations between the model and the 
ab-initio results, both data sets clearly follow the same 
trend. This result further supports TVS as a spectro- 
scopic tool, but it also underlines that quantitative in- 
formation about the molecular level position can only be 
obtained from TVS provided some knowledge of the junc- 
tion asymmetry. The calculated numbers for the data 
points are listed in Fig. [3l 

The current at finite bias voltage is calculated us- 



ing DFT in combination with a non-equilibrium Green 
function (NEGF) method, fohowing the principles of 
Ref. [IBj. Our DFT-NEGF method is implemented in 
GPAW, which is a real space electronic structure code 
based on the projector augmented wave method lid llTj. 
We use the PBE exchange-correlation functional [18| , and 
a 4 X 4 /c-point sampling in the surface plane. The elec- 
tronic wave functions are expanded in an atomic orbital 
basis [17|. In all calculations, the molecule and the closest 
Au layers are described by a double-zeta plus polariza- 
tion (dzp) basis set, while the remaining Au atoms are 
described by a single-zeta plus polarization (szp) basis. 

The upper panel of Fig. [3] shows the atomic structure 
of two representative junctions: a hexanethiol (left) in 
an STM-like configuration with the STM-tip ~ 4 A away 
from the molecule, and anthracenethiol (right) between 
two fiat Au (111) surfaces with a distance of 1.5 A be- 
tween the right Au electrode and the closest H atom 
of the molecule. For all STM-setups, we initially re- 
lax the molecule on a single Au(lll) surface and sub- 
sequently add the STM electrode without further relax- 
ations. For the systems with two flat Au surfaces, we 
relax the molecule and the two closest Au layers. Ex- 
cept for the benzene junctions (S-Ben-S and HS-Ben- 
SH), the linking S atom relaxes into a bridge site of the 
Au(lll) surface shifted slightly towards the hollow site. 
The structure of the symmetric benzene-dithiol (S-Ben- 
S) junction is taken from Ref. [19| where the molecule 
binds to an Au adatom together with a SCH3 unit. In the 
other benzene junction (HS-Ben-SH) the hydrogenated 
sulfur atoms bind to Au adatoms of the Au(lll) surfaces 
following Ref. [2o|. 

In Fig. [3] we also list for each junction the obtained 
TVS minimum voltage, Fmin, the position of the HOMO 
level relative to Ep^ and the calculated asymmetry pa- 
rameter. Figure |4] illustrates how we determine these 
quantities, in the case of the C6-S junction. The TVS 
minimum is simply found from the minimum in the 
Fowler- Nor dheim plot, as show in the inset. To deter- 
mine the HOMO level, we project the density of states 
onto the carbon atoms. Molecular levels then appears 
as clear peaks in the projected density of state (PDOS) 
vs. energy plot as shown in Fig. HI We note that there 
are situations where a peak in the PDOS does not show 
in the transmission function due to symmetry mismatch 
between the molecular orbital and the electrode wave- 
functions. This is indeed the situation for Ph-S, Naph- 
S and Anth-S in the STM setup where a peak in the 
PDOS is observed at E — Ep ~ -0.55eV, but only a 
vanishingly small transmission peak occurs at this energy. 
In these three cases we take the second PDOS peak at 
E - Ef = -2.35eV as the HOMO level. At this energy, 
the Au tip wavefunctions also have d-character, which 
enables non-zero coupling to the molecular 7r-orbital. 

At finite bias voltage, the energy of the HOMO level 
is shifted by the bias voltage, (c.f. Fig. [1]). Figure 
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Figure 3: Top: Atomic structure of two of the major groups 
of molecular transport junctions investigated in this work and 
in experiments 2_, ^, _5] , namely alkanes in an STM-like con- 
figuration (top left), and anthracenethiol (Anth) between two 
flat Au (111) surfaces (top right). The table in the lower part 
shows schematically all the considered molecular junctions to- 
gether with calculated values of Vmin, sh — Ep, and 77. The 
electrode configurations are indicated with gray boxes show- 
ing both the flat Au surface (rectangles) and STM-like geome- 
tries (pentagons). The super-script ^ indicates that the data 
are obtained from non-self consistent calculations by manually 
setting 7] — 1/2. The HOMO energy, sh, refers to the first 
peak in the Kohn-Sham (PBE) projected density of states. 
See text for more details. 



m shows the PDOS of the Ce-alkane junction under bias 
voltage y = V (red solid) and V = 2.0 V (blue dashed). 
The equilibrium HOMO energy of en — Ep = — 3.6eV 
shifts upward by A = 0.4 eV when the bias is applied. 
The asymmetry factor, 77, is calculated from the shift, A, 
according to 77 = A/V = 0.2. 

Notice that even for the very asymmetric STM config- 
urations, the asymmetry factor is never larger than 0.35. 
This is due to the relatively small distance between the 
molecule and the Au tip (4.0 A). Increasing this distance 
would increase r] further, however, this is computation- 
ally problematic due to the finite range of the atomic or- 
bital basis. Instead we simulate a completely asymmetric 
junction (77 = 0.5), by using the zero-bias transmission 
function and fixing the electrode chemical potentials to 
fiL = Ep and (ir = Ep — V. The data for the three last 
junctions in Fig. [3] are obtained in this way and are thus 
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Figure 4: Projected density of states on the C-atoms in the 
Ce-alkane junction shown in Fig. [S] The sohd red hne 
refers to Vbias = OV and the dashed blue hne refers to 
Vbias = 2.0 V. The HOMO level (at Vbias = OV) is located 
dX E — Ef — — 3.6eV. We determine the asymmetry factor, 
?7, from the shift, A, of the HOMO level with the bias volt- 
age. The inset shows the calculated Fowler- Nor dheim plot 
with the characteristic TVS-minimum at Vmin = 2.4 V. 



non-self consistent results. 

The deviations of the numerical data from the ana- 
lytical result in Fig. [2] may be due to several reasons, 
(i): The ab-initio transmission functions do not have a 
perfect Lorentzian shape as assumed in the analytical 
model, (ii) The dependence of the HOMO energy on the 
bias voltage is not strictly linear, (iii) The bias voltage 
influences the shape of the molecular orbitals which in 
turn affects the coupling strength, (iv) It may in some 
cases be difficult to accurately determine £h and r^, as 
discussed above. 

The TVS minimum voltage was originally interpreted 
as a transition from tunneling to field emission [J] . How- 
ever, according to our calculations, Vmin does not mark 
such a transition, and in line with Refs. [ll|, |2l| we do 
not view Vmin as a special transition voltage. Rather the 



usefulness of TVS relies on the direct relation between 
Fmin and Eh when the transmission function can be well 
described by a Lorentzian. 

In the conducting AFM measurements of Ref. 4, the ra- 
tio \sH — EF\/Vmin was found to vary between 2.1 and 2.6 
for a set of molecules which included Naph-S, BP-S, and 
TP-S also considered in this work. This is significantly 
larger than the value of 1 predicted by Simmons bar- 
rier model. According to our model the ratio should be 
smaller than 2.0 which is obtained for a completely asym- 
metric junction ij] = 0.5). To resolve this puzzle, we first 
note that our ab-initio calculations yield r] ~ 0.3 for the 
three molecules (c.f. Fig. [3]) giving \eH — EF\/Vmin ^ 1-6 
according to our model. Next, we note that the pres- 
ence of the AFM tip in the transport measurements will 
lead to a renormalization of the HOMO energy due to 
an image charge effect [22[ which is not present in the 
UPS measurements. Based on the method described in 



Ref. |23| we estimate the image charge interaction to be 
0.5 — 0.35 eV depending on the size of the molecule. Cor- 
recting the UPS values for en by these values lead to 
ratios |^i^ — £^F|/Knin in the range 1.6 — 1.8 in good agree- 
ment with the model prediction of '^ 1.6. 

Finally, we note that the well known inability of DFT 
to describe energy gaps and level alignment of molecules 
at surfaces |22| does not affect the conclusions of the 
present work. This is because, according to Fig. [TJ the 
ratio {eh — £^F)/^in is independent of the value of €h- 
In particular, the dependence of {eH—EF)/Vmin on asym- 
metry is expected to be a general result independent of 
the absolute position of the molecular levels. 

In conclusion, we have performed extensive ab initio 
DFT transport calculations to simulate transition voltage 
spectroscopy for a large number of molecular junctions. 
The numerical data closely follow the trend expected 
from an analytical model with a Lorentzian shaped trans- 
mission function. We have explicitly shown that in order 
to use TVS as a quantitative spectroscopic tool to probe 
the molecular levels, it is necessary to take the asymme- 
try of the molecular junction into account. The present 
analysis should therefore be considered in future applica- 
tions of transition voltage spectroscopy. 
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